Enhanced heat transport by turbulent two-phase Rayleigh-Benard convection 
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We report measurements of turbulent heat-transport in samples of ethane (C2H6) heated from 
below while the applied temperature difference AT straddled the liquid-vapor co-existance curve 
T^{P). When the sample top temperature Tt decreased below T^, droplet condensation occurred and 
the latent heat of vaporization H provided an additional heat-transport mechanism. The effective 
conductivity Ae// increased linearly with decreasing Tt, and reached a maximum value A*jj that 
was an order of magnitude larger than the single-phase Ae//. As P approached the critical pressure, 
A*jj increased dramatically even though H vanished. We attribute this phenomenon to an enhanced 
droplet-nucleation rate as the critical point is approached. 

PACS numbers: 47.27.te, 47.20.Bp, 47.55. pb, 47.55. db 



Turbulent heat transport in a fluid heated from below 
(Rayleigh-Benard convection or RBC) has been a topic of 
intense fundamental research for some time [l|, Q . Most 
of this work was done under conditions where the fluid 
was in a single phase, that is far from any thermodynamic 
phase-transitions. Here we focus on the case where the 
applied temperature difference AT = Tb — Tt (Tb and Tt 
are the bottom and top temperatures respectively) spans 
a liquid- vapor phase-transition line Tif,{P). In that case 
condensation and vaporization ("boiling") can provide 
heat-transport mechanisms in addition to the usual tur- 
bulently advected heat current. For these mechanisms 
the latent heat of vaporization H clearly plays a ma- 
jor role. This case is of great practical importance be- 
cause the exceptionally high effective thermal conductiv- 
ity Ae// enables numerous industrial applications, includ- 
ing for instance miniaturized heat exchangers and per- 
formance enhancement in process industry, ft has been 
studied extensively from an engineering viewpoint, and 
a wealth of empirical correlations based on these studies 
has been used in designs that range from miniaturized de- 
vices for coolirig of computer components to large-scale 
power plants. [1, 0, 0| 

Here we report on a study that was intended to ad- 
dress some of the fundamental physical aspects of this 
problem. We measured Ae/ / in cylindrical samples of di- 
ameters D about equal to their heights L using ethane 
(C2H6) below but near its critical point (CP). The sample 
pressure P and AT were held constant while the mean 
temperature Tm = {Tt+Tb)/2 was changed in steps much 
smaller than AT through the two-phase region. We fo- 
cused primarily on the parameter range where the heat- 
transport enhancement SXeff was due to condensation 
near the top plate where the system had a large thermal 
gradient and where the top temperature was below T^ 
while most of the sample remained in the vapor phase at 
temperatures above T^. As Tt was gradually lowered be- 
low T^, Ae// increased linearly above the single-phase 
value. In this regime Ae// was completely reproducible 



and independent of history. Shadowgraph images showed 
that condensation was by droplet formation rather than 
by film condensation. 0] Even though H vanishes at 
the CP, the largest enhancement 5A*yy (P) of SXeff [P, T) 
at a given pressure increased rather dramatically as P 
approached the critical pressure Pcp- Since the heat 
transport depends not only on H but also on the rate of 
droplet formation, this result implies a droplet-nucleation 
rate that increased more rapidly than H decreased as 
P — > Pcp- This result is qualitatively consistent with 
classical nucleation theory [7|; but clearly our system, 
with a large thermal gradient just below the top plate 
and with vigorous fluctuations, is more complicated than 
those treated before by that or more advanced theories. 
Interestingly, SXeff was essentially the same for two dif- 
ferent samples, one with a finely machined copper top 
plate and the other with an optically flat sapphire top 
plate, suggesting that surface roughness did not influence 
the nucleation rate signiflcantly and that the nucleation 
process was homogeneous. 

The apparatus had been used for several previous in- 
vestigations of turbulent RBC 0, i, [H El, [13. Se- 
quentially we installed two different high-pressure sample 
cells. The flrst, cell A, was a cylinder with T=7.62 cm 
and D=7.63cm. It had been used for turbulent heat- 
transport measurements in gases [13, [l3|- Its top and 
bottom consisted of thick copper plates, with finely ma- 
chined inner surfaces, that fit closely into a side wall 
made of high-tensile-strength stainless steel. For flow- 
visualization with the shadowgraph method [l5| we used 
cell B which had D = 10.16 cm and L = 9.84 cm. 
It had an optically flat sapphire top plate but a cop- 
per bottom plate with a polished surface and an evap- 
orated gold film that served as a mirror. The measure- 
ments of Ae// reported here are for cell A, but cell B 
yielded largely equivalent results. Heat was applied at 
the sample bottom by a metal-film heater covering the 
bottom-plate area uniformly. The top plate was cooled 
by a circulating water bath. Both Tf, and Tt were held 
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FIG. 1: Pressure-temperature plane showing the coexistence 
curve T^{P) (solid line) and the CP (circle) of C2H6. Dashed 
lines: some of the isobars used here. 
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FIG. 2: Effective conductivity Xeff as a function of the mean 
fluid temperature Tm (lower abscissa) and the temperature 
fraction = 0.5 + {Tm - T^{P))/AT (upper abscissa) for 
P = 44.77 bars and AT = 1.00 K. Triangles: Xeff is time- 
independent. Circles: Xeff is time-dependent (the time av- 
eraged values are shown; typical fluctuation amplitudes are 
indicated by the bar at Tm ~ 28.03° C). Solid (open) sym- 
bols: Tm was increased (decreased). 



constant within a milli-Kelvin or better. The sample 
was connected to an external volume through a capil- 
lary. The temperature of this volume was controlled in 
a feedback loop with a pressure gage so as to hold the 
sample pressure constant within 10~^ bars. C2H6 has 
well known properties 17, 1^ and a conveniently located 
CP at Tcp = 32.172°C, Pcp = 48.72 bars. Its phase- 
separation curve T^{P) is shown in Fig.[TJ 

In Fig. [2] measurements of Ae// = QL/AT {Q is 
the heat-current density) at constant P = 44.77 bars 
and AT — 1.00 K are plotted as a function of T„i 
(lower abscissa) or of the temperature fraction = 
0.5 (T™ - T^{P))/AT (upper abscissa). For this pres- 
sure T^{P) = 28.148° C. The two vertical dashed lines 
show the temperatures r„ = T^{P) ± Ar/2. For (/> > 1 
the entire sample was in the vapor phase and Ae// was 
close to values obtained from an extrapolation of mea- 
surements made before 13|,|lJ| but further away from T^. 
For (j) < the entire sample was in the liquid phase, and 
again Ae// was consistent with other measurements. As 
(j) was lowered from (f> > 1 into the two-phase region below 
(j) = 1, Xeff initially increased linearly as a function of (j) 




FIG. 3: (a): The shift T^ - Tt,on of the onset of heat- 
current enhancement as a function of the applied tempera- 
ture difference AT. The solid line has a slope of 0.043. (b): 
Xeff as a function of <f). Solid circles: AT = 2.00K. Open 
squares: AT = 5.00K. Triangles: AT = 8K. Solid squares: 
AT = lO.OOK. Open diamonds: AT = 15.00K. The pressure 
was 43.04 bars. 



and reached a maximum value A*jy at (jf ~ 0.43 that was 
nearly an order of magnitude larger than in the vapor. 
There was a sharp onset of this heat-transport enhance- 
ment, but there was no discontinuity. The solid (open) 
symbols correspond to data taken with increasing (de- 
creasing) Tm or 4>. One sees that Ae/ / was independent of 
this past history and highly reproducible. This region of 
linear increase, we shall show, corresponded to a sample 
filled with vapor but with droplet condensation occurring 
at the top plate. As (j) dropped slightly below 0* , the heat 
transport became time dependent and the time-averaged 
values became history dependent. Further reduction of (p 
(but still with (j) > Q) led once more to time-independent 
states, but with relatively small heat-transfer enhance- 
ments that also varied from run to run. This parameter 
range corresponded to a sample filled with liquid and 
with vaporization (or "boiling" ) occurring at the bottom 
plate[16l|. The transition from a mostly vapor-filled to 
a mostly liquid-filled state near (/) = (/)* could easily be 
seen in the experiment because it led to a discontinuous 
increase of the temperature of the external volume used 
for the pressure regulation. 

An interesting aspect of the onset of the heat-transport 
enhancement is that it occurred at Tt,on < T^{P). The 
shift T^ — Tt.on increased roughly linearly with AT, as 
shown in Fig. [3^ for a pressure of 43.04 bars. This linear 
increase implies a constant shift of the temperature frac- 
tion 4>on at onset below = 1 ; it is illustrated in Fig. [3)3 
where 1 — (f)on — 0.05 independent of AT. 

In Fig. UK one sees that a larger AT leads to a slower 
increase of Ae// with decreasing Tt. Figure SJj reveals 
that, for sufficiently large T^ — Tt, the data for the heat- 



FIG. 4; (a): The effective conductivity Ae// and (b): Q at 
a pressure of 43.04 bars and at several applied temperature 
differences AT. Open circles: AT = LOOK. Solid circles: 
AT = 2.00K. Open squares: AT = 5.00K. Solid squares: 
AT = lO.OOK. Open diamonds: AT = 15.00K. Solid line in 
(b): Q = -0.212(rt - T^) - 0.013(T - T^^f W/cm^ 



current density Q approach a single curve, showing that 
the heat transport is determined by — Tt and not by 
AT. 

In Fig. [5] we show shadowgraph images taken in cell 
B. It is difficult to interpret these images quantitatively 
because they give a vertical average and the vertical lo- 
cation of any feature remains unresolved. Nonetheless 
they give useful qualitative information. Image (a) for 
(j) = 1.01 is for the single-phase vapor region. At this 
point the Rayleigh number Ra = (3g/^TL'^ / [nv] {g is the 
gravitational acceleration and /9, and k are the ther- 
mal expansion coefficient, kinematic viscosity, and ther- 
mal diffusivity respectively) was about 3.8 X 10^° and the 
Nusselt number Nu — Xe/f/X (A is the diffusive thermal 
conductivity) was close to 220. The structure seen in 
the shadowgraph corresponds to plume activity and fluc- 
tuations of the highh^ turbulent single-phase system. A 
movie for this case [19|] shows the turbulent time depen- 
dence and reveals the existence of a large-scale circula- 
tion which swept the structures along [20|, [21 1. Image 
(b), for (j) = 0.92, is in the two-phase region where a 



thin layer of fluid in the top thermal boundary layer had 
been rendered meta-stable. One sees a new feature: there 
were numerous small dark circles which we interpret to 
be liquid droplets. They moved laterally and were swept 
toward the side wall, we believe by the prevailing large- 
scale circulation. A movie of this case can also be found 
elsewhere [3]. Images (c) and (d) are for (j) — 0.85 and 
0.71 respectively, and the abundance of droplets is seen 
to have increased as (f> decreased. 

Images (e) and (f ) are for < (p < (p* where the sample 
is mostly liquid-filled. Somewhat above cf) = boiling 
started but Ae// was not enhanced very much. A small 
number of isolated gas bubbles can be seen in (c), with 
one of them identified by the small white arrow in the 
lower left part. The bubbles meander chaotically in the 




FIG. 5: Shadowgraph images from cell B at different mean 
temperatures with (a) (j) = 1.01, (b) 0.92, (c) 0.85, (d) 0.71, 
(e) 0.32, and (f) 0.58. For this experiment P = 41.37 bars, 
AT = 0.50/C, and Tm was increasing. Movies for (j> — 1.01 
0.92, 0.32, and 0.58 are available elsewhere [131. 



lateral direction and, upon rising, re-dissolve in the cooler 
sample interior. At larger cf> < (j)* , the small bubbles 
collect in one large bubble located under the top plate as 
shown in (f). Inside that bubble condensation is taking 
place as evident from the many drops that form within 
it. These processes are illustrated better by the movies 
[lit for (5e) and (5f). In the boiling range, with (jj < (jj* , 
measurements of Ae// were irreproducible from one run 
to another. 

The latent heat vanishes at the CP, and above the 
CP the liquid and gas states become indistinguishable. 
Thus one might expect initially that the heat-transport 
enhancement in the two-phase region should be weakened 
as the CP is approached and cease to exist for P > Pcp- 
In Fig. [5] we show measurements of Ae// with AT = O.IK 
at different P. Remarkably, as P — > PcPj Kiff increased 
dramatically. As P exceeded Pcp, A*yy decreased again. 

At the CP (3 diverges and k vanishes. As a result 
Ra and the Prandtl number Pr = v/n become infi- 
nite. For these rapidly varying conditions we estimated 
the heat transport contributed by single-phase turbulent 
Rayleigh-Benard convection Xrb using the predictions 
from Ref. [22] with fluid properties evaluated at Tm- 
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FIG. 6: The effective conductivity Ae//, measured with AT = 
0.1 K, as a function of Tm for different isobars. From left to 
right, the data are for P=46.14, 46.83,47.52, 48.21, and 49.25 
bars. The solid curves are estimates for single-phase turbulent 
convection using the properties at Tm and the predictions of 
Ref. [2^ . The dotted line shows the critical temperature. 
The inset shows the maximum excess contribution SX^jf = 
A*^ f — \*iiB from the nucleation process to Ae/ / as a function 
of the reduced pressure e^. 

These estimates are shown by the sohd curves in Fig. [51 
For P above Pcp there is reasonable agreement between 
the prediction and the measurements. Below Pcp, as 
expected, the estimate is much too small and there is 
a much larger contribution from condensation or boiling. 
The inset of Fig.lHlshows the excess 5\*fj = i^lff-^RB) 
as a function of the reduced pressure Cp = (Pcp—P) /Pcp 
i^*R.B maximum value of XpB at a given P). The 

maximum heat transport contributed by the nucleation 
process increased by an order of magnitude as the fluid 
pressure approached Pcp- 

In this Letter we reported on heat-transport mea- 
surements by turbulent Rayleigh-Benard convection of 
ethane under conditions where the applied tempera- 
ture difference AT straddled the liquid- vapor coexistence 
curve T^{P). As the top temperature Tt was lowered 
quasi-statically below T^, the effective conductivity Ae// 
was enhanced by droplet condensation at the sample top. 
The droplet formation was observed by shadowgraphy. 
With decreasing Tt and starting at Tt,on < T^, X^ff ini- 
tially increased linearly and reached a maximum at 
Tj* that was an order of magnitude or more larger than 
typical values in the single-phase regions. Here Ae// was 
reproducible and history independent. The shift of the 
onset T^ — Tt^on increased roughly linearly with AT. 

The maximum enhancement <5A*y^ above the single- 
phase value increased dramatically as the pressure ap- 
proached the critical value Pcp- Since the latent heat 
H vanishes at the critical point, this implies that the 
droplet-nucleation rate increased at a sufficiently large 
rate to more than overcome the diminished contribu- 



tion from H. We do not know of a droplet-nucleation 
theory that would be applicable quantitatively in the 
presence of the steep thermal gradient and the vigorous 
fluctuations characteristic of turbulent convection. How- 
ever, under the more benign circumstances of an isother- 
mal meta-stable fluid classical nucleation theory sug- 
gests that the nucleation rate should be proportional to 
exp(— AF/ZcsT), and that the difference in free energy 
AF between the vapor and the droplet should vanish at 
the critical point. Qualitatively this implies an enhanced 
nucleation rate as P ^ Pcp- 
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